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Resonant-pulse operations on the buried 
donor charge qubits in semiconductors 

L. A. Openov 

Moscow Engineering Physics Institute (State Univer- 
sity), 115409 Moscow, Russia 

A new scheme is proposed for rotations 
of a double-donor charge qubit whose logical 
states are defined by the two lowest energy 
states of a single electron localized around 
one or another donor. It is shown that mak- 
ing use of the microwave pulses tuned to the 
resonance with an auxiliary excited molec- 
ular level allows for implementation of vari- 
ous one-qubit operations in very short times. 
Dccoherence effects are analyzed by the ex- 
ample of the P^:Si system and shown to be 
weak enough for experimental realization of 
this scheme being possible. 



the logical states are formed by the lowest two energy 
states of the remaining valence electron localized at the 
left or the right donor, |0) = \L) and |1) = \R), see Fig. 
1. Initialization and readout of the qubit are facilitated 
by a single electron transistor. The surface electrodes are 
used to control the qubit through the adiabatic variations 
of the donor potentials. 
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FIG. 1: The logical states |0> 
buried donor charge qubit. 
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PACS Numbers: 85.35.-p, 03.67.Lx, 73.20.Hb 



Semiconductor-based devices seem to be very promis- 
ing for a scalable quantum computing technology 0. 
The qubits can be encoded, e. g., in the nuclear or 
electron spin states 0, HI • Although the relatively long 
coherence times make the spin-based qubits good can- 
didates for quantum computation, the single-spin mea- 
surement still remains a significant challenge [3J, |j] . The 
charged-based qubits in semiconductors are currently dis- 
cussed as well, their logical states being encoded in the 
orbital degrees of freedom of an electron occupying the 
quantum-dot structure [M IE 01 HI • Coherent oscillations 
of the double-dot qubit have been observed Q, and the 
readout schemes have been proposed [E]. In spite of the 
fact that decoherence of the charge-based qubits [lfj is 
much stronger than of their spin-based counterparts, the 
charge-based qubits are nevertheless believed to be real- 
izable at the present technological level due to their short 
operation times. 

In view of the technical difficulties concerning manu- 
facturing of the quantum dots with predetermined char- 
acteristics, it may appear more reasonable to make use 
of natural atoms (instead of " artificial" ones) as the lo- 
calization centers for the electrons carrying qubits. Re- 
cent advances in atomically precise placement of single 
dopants in semiconductors 0, 0] make possible the 
construction of solid-state atomic qubits. Hollenberg et 
al. proposed a two-atom scheme where the charge qubit 
consists of two dopant atoms beneath the semiconductor 
surface [l3.lF3l |. One of the donors is singly ionized, and 



Here we propose an alternative scheme for rotations 
of the buried donor charge qubit. Instead of applying 
biases to the surface gates, the qubit is driven by two 
microwave pulses tuned to the resonances between the 
localized states \L), \R) and one of the excited molecular 
states delocalized over the double-donor structure. We 
present the analytical solution for the three-level model of 
the unitary electron evolution and show that, depending 
on the specific values of frequencies, phases, amplitudes, 
and durations of the pulses, various one-qubit operations 
can be implemented in times as short as ~ 1 ps, i. e., 
more than one order of magnitude shorter than in the 
original proposal 0,0. At such times, decoherence is 
sufficiently weak, making it possible at least to investi- 
gate the small-scale devices and thus to demonstrate the 
experimental feasibility of the scheme. 

We consider a singly ionized pair of dopant atoms in a 
semiconductor host by the example of the P^:Si system. 
The Hamiltonian of the remaining valence electron is 



H Q = J2 E n\Xn)(Xn\ 



(1) 



n=l 



where E n and \xn) are, respectively, the one-electron 
eigenenergies and eigenstates of the diatomic ion. In 
order to avoid the sophisticated numerical calculations 
|l4j , we ignore the conduction-band anisotropy, the inter- 
valley terms, the surface effects, etc., and make use of a 
simple effective mass approximation that is comm only 
used for the semi-quantitative considerations 0, llfij . 
Then, in the case that the two donors are equivalent, 
the problem reduces to that for the hydrogen-like molec- 
ular ion with the effective Bohr radius a* R a: 3 nm and 



2 



the effective Hartree unit of energy E* = e 2 /ea* I 



40 



meV, where e ~ 12 is the dielectric constant for sili- 
con. The lowest states and |% 2 ) are the molecular 
states lsug and 2pa u , whose wave functions are, respec- 
tively, symmetric and antisymmetric about the midpoint 
of the line joining the two donors. At large donor sepa- 
ration, R >> a* B , one has E\ w E 2 « —E*/2, while the 
difference A £21 = £2 — £1 is exponentially small [l7j . 
AE21/E* = Axe-*- 1 [1 + 0(1/2;)], where x = i?/a^. 
The excited state \xa) with the energy £3 « — £*/8 is 
well separated from the states |xi) and \x%) by AE31 « 
A£ 32 w 3£*/8 « 15 meV. At J?/a| > 6 this is the 
molecular state 3rfcr g . If the qubit is biased by the static 
gate voltages, so that the difference Er — El in the ener- 
gies of the lowest states \L) and \R) localized at different 
donors greatly exceeds the value of AE21, then the states 
\L) and \R) are good approximations to the lowest eigen- 
statcs, and 



H = E L \L){L\+E R \R){R\+J2 E n\Xn){Xn 



(2) 
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Let the qubit interact with an external electromagnetic 
field E(t) = Eoi(t) cos(uj L t) + F, 02 (t) cos(uj R t + <p) that 
has two components oscillating at frequencies ojl = (E3 — 
E L )/h and u) R = (£3 - E R )/h, where E i(i) and E02W 
are the slowly varying envelopes. Then the Hamiltonian 
becomes 



H(t) = H a + V(t) . 
The interaction term V(t) is 



V(t) = E(t) 



d L \ X 3)(L\ + d R \ Xa )(R\ + h.c. 



(3) 



(4) 



where = (xa| — er|L, i?) are the electric dipole mo- 
ments for the transitions \L,R) ^ |X3)- I n the resonant 
approximation one has 



V(t) 



+ -e--«*-^A fl (i)|x3><i?|+/i.c. 



(5) 



where A L (i) = E i(t)d L and Ar(*) = E 2(*)dfl. Here 
we restrict ourselves to the rectangular pulse shapes, so 
that both E01 (t) and E02 (t) are constant at < t < r op 
and zero elsewhere. 

It is straightforward to solve the non-stationary 
Schrodinger equation 



Ot 



(6) 



for the state vector 



|*(t)) = C L {t)e- lELt / h \L) + C R {t)e- %ERt ' h \R) 
+C 3 (t)e-* E ^ h \ X3 ) (7) 



and to find the coefficients Cx(t), C R (t), and C3(i) pro- 
vided that |*(0)) = a\L) + f3\R): 



C L (t) = a 

-a 



1 - 



|Al| 2 + M 2 

2\l\ R e-^ . 2 



sin 2 (fit) 



A L | 2 + |A fi 



2 sin 2 (m) , 
2AiA| ? e 



'|Al| 2 + |a«| 



sin 2 (fti) 



+/3 
<2>(t) 



1 - 



2|A*| : 



|AL| 2 + |A fl | : 



■sin 2 (fit) 



.aA L + 0X R e 



-i(f> 



IAhI 2 



sin(2fii) 



(8) 



where f2 = \J\X L \ 2 + \X R \ 2 /Ah. From Eq. © one can 
see that at t = r op — nk/20, (k is a positive integer) the 
coefficient C3 vanishes, and the state vector |^(t)) re- 
mains in the qubit subspace {\L}, \R}}. So, the auxiliary 
excited state |x3) assists the qubit evolution by means of 
the electron transfer between the states \L) and \R) as 
the drivin g fie ld is on but remains unpopulated after the 
field is off ]3. 

It follows from Eqs. J7J and JSJ that the relative phase 
shift operation, 



— iE L T op /h 



\^(rop)} = e 



is implemented at r op 
eration, 



t\L)+pt 



-i(E R -E L ) Top /h 



\R) 



(9) 

irk/Q. The quantum NOT op- 



*(r op )) = ±e 



-iEl, T p/h—i(j> 



0\L)+a\R) 



(10) 



is realized at r op — ir(2k — l)/2f2 if Xl = T^R an d 4> — 
Ttn + (E R — EL)T op /2h (n is an integer). The Hadamard 
transformation, 



,/ft 



a + (3 



a-0 

~7f 



\R) 



, (11) 



is performed at r op — 7r(2fc — l)/2fi if [E R ~ EL)r op /h — 
2irm (m is a positive integer). The plus sign in Eq. (jl 1|> 
corresponds to <p — 2wn and Xl = — X R (y2 — 1) or <fi = 
n(2n + 1) and Xl = X R (y2 — 1), and the minus sign 
corresponds to <fi = 2irn and Xl = X R (\/2 + 1) or cj> — 
■n(2n + 1) and X L = -X R (y/2 + !)■ 

For the field amplitudes Eq ~ 1 V/cm, the operation 
time is r op ~ I/O ~ ft/IAi^J ~ h/ea* B Eo ~ 1 ns. In- 
crease in the pulse intensity will cause the value of r op 
to decrease down to the picosecond time scale, so that 
the value of T op can be made much shorter than in the 
case that the qubit is manipulated by adiabatically vary- 
ing the potentials of the surface gates Note that 
the energies El and E R should be sufficiently different 
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from each other, Er — El ~ 1 meV, in order the qubit 
rotations could be implemented in times r op ~ 1 ps. 

The uncontrolled interaction of the quantum system 
with its environment leads to entanglement between the 
states of the system and the environmental degrees of 
freedom. This disturbs the unitary evolution of the sys- 
tem and results in the loss of coherence. There are var- 
ious sources of decoherence in solids. For the buried 
donor charge qubit decoherence due to the phonon emis- 
sion/absorption processes was studied in Refs. 0, 
and found to be much weaker than decoherence due to 
both Nyquist- Johnson voltage fluctuations in the surface 
electrodes and 1 /f noise from the background charge fluc- 
tuations. Contrary to this statement, here we show that 
phonons are the main cause for decoherence at short op- 
eration times. For simplicity, we consider the qubit at 
zero temperature and assume isotropic acoustic phonons 
with the linear dispersion law, w q = sq, where s is the 
speed of sound. 

Electron-phonon coupling in confined systems is de- 
scribed by the Hamiltonian [2£j 



H ep = A (q)p(q) 



b+ + b_ 



(12) 



where fo+ and 6 q are, respectively, the operators 
of creation and annihilation of a phonon with the 
wave vector q, p(q) — J <ire iqr /3(r) is the Fourier 
transform of the electron density operator p(r) = 

Smn v f / m( r )*n( r )l m )( n lj an d Hq) ^ s tne microscopic 
electron-phonon interaction matrix element that can be 
expressed in terms of the deformation potential D and 

1/2 

the density of the crystal p as A(q) = qD (h/2pLo o p) 1 , 
with fl being the normalizing volume. 

Since at R >> a* B the overlap (L\R) between the 
orbitals (r|L,J2) = {n(a B ) 3 )-^ 2 exp(-|r - T LjR \/a* B ), 
where vl.r. = ^f(R/2)e x are the donor coordinates, is 
negligibly small, the transitions \L) ^ \R) are sup- 
pressed, and decoherence of the lowest states \L, R) is 
entirely due to dephasing processes, so that the diagonal 
elements of the density matrix remain unchanged, while 
the off-diagonal elements are 0] 



PLB{t) = p LR (0)e- B2 ^ E «- E ^ 



(13) 



where 



B2(t)= 8 ^| 5 (q)| 2 ,, iL ,,. q 



h 2 ^ ujI 



U! n t 



, (14) 



and 



s(q) 



A(q) 



-»A(q)- 



[{L\e l ^\L) - (i?|e lqr |i?) 
sm(q x R/2) 



l + {qa B f/A\ 



(15) 



At T op > a* B / s one has 

B 2 (T op ) : 



D 2 



3ir 2 phs 3 (a* B ) 2 ' 



(16) 



so that the spectral function l|14|) appears to be a material 
constant, being about 6 • 10 -3 for P^~:Si, and the error 
rate (i. e., the error generated during the operation time) 
is 



D{t op ) 



3-10" 



(17) 



Since the excited level |x3) becomes temporarily pop- 
ulated during the resonant-pulse operations on the P^~:Si 
qubit, the phonon emission processes \\3) — ► \L,R) also 
contribute to decoherence at T = 0. For the dou- 
ble donor orientation along the x-axis one has w 
[\2S) L ~ \2P x )l + \2S) R + \2P x ) R ]/2 at R » a* B 
and Er — El « E3 — El,r- Neglecting the expo- 
nentially small overlap between the localized atomic- 
like orbitals (t\L,R), (t\2S) l ,r = (87r(a/k) 3 ) _1/2 (l - 
|r - T L>R \/2a B ) exp(-|r - v L ,R\/2a B ), and (r\2P x } L , R = 
(32TT(a B ) 5 )-^ 2 (x-x L ,R)ex V (-\i 



2a* B ) centered 



at different donors, we have 



( X3 |e lqr |L,i?> = 2^. 



(qa* B ) 2 T *i(?xOfl) 



(qa B ) 2 } 



(18) 



so that the relaxation rate at T = is 12JJ 



T =T £ \K<l)\ 2 \{Xs\e^ T \L,R)\ 2 5(hLJ -nw q 

q,L,R 



8D< 



Trphs 2 (a* B 



* \5 4 



* + (Qoa* B ) 2 



7* ^2 



(19) 



where hu> a = hq s = E 3 — E L ~ E 3 — Er, ps 3E* /8. From 
Eq. one has T » 3 • 10 7 s" 1 for P+:Si. We see that 
at T op < 100 ps the error rate due to the phonon emission 

0, 



processes 



D(t op ) = 1 



,-rr 



(20) 



is lower than the value of D(r op ) due to dephasing. So, 
the phonon-induced error rate at T = and short oper- 
ation times is D(r ov ) m 3 ■ 10~ 3 . At finite temperatures, 
such that k B T > hwo, where Hloq = hs/a* B « 2 meV for 
dephasing processes and fiwo = E 4 — E 3 « Er — El ~ 1 
meV for the processes of the phonon absorption by an 
electron temporarily occupying the excited state, the er- 
ror rate increases by a factor of ~ k B T/hujQ, i. e., changes 
slightly at T < 10 K. 

The error rate due to phonons should be compared to 
the error rates due to other sources of decoherence. The 
lowest bounds for the decoherence times associated with 
the Johnson noise from the gate s and the environmental 
charge fluctuations are [IS Il3l ll5T | r ~ 1 [is and r ~ 
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1 ns, respectively, so that the corresponding error rates 
D(t op ) = 1 — exp(— t op /t) do not exceed that due to 
phonons at r op < (1 10) ps. Hence, the performance 
of the buried donor charge qubit appears to be limited 
primarily by the electron-phonon interaction. 

In conclusion, we proposed a scheme for fast rotations 
of the buried donor charge qubit through an auxiliary- 
state-assisted electron evolution under the influence of 
the resonant microwave pulses. This scheme allows for 
implementation of the one-qubit operations in times as 
short as r op ~ 1 ps. By the example of the P^Si qubit, 
we have shown that dephasing due to acoustic phonons 
is the main source of decoherence at T < 10 K and op- 
eration times T op = (1 -j- 10) ps. The error rate is about 
3 ■ 10 -3 , i. e., greater than the fault-tolerance threshold 
for quantum computation |22| but low enough for co- 
herent qubit manipulation being possible, at least in the 
proof-of-principle experiments on one-qubit devices. The 
coupling of the double-donor qubits via the Coulomb in- 
teraction allows, in principle, to realize the conditional 
gates. It would be also worthwhile to search for other 
materials and/or doping elements for the buried donor 
charge qubits, in order to weaken the decoherence effects. 
Although we restricted ourselves to rectangular shapes of 
the resonant pulses, our consideration can be generalized 
to other pulse shapes 

Discussions with A. V. Tsukanov, L. Fcdichkin, and 
M. S. Litsarev are gratefully acknowledged. 
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